Revised WMAP constraints on neutrino masses and other extensions of the minimal 

ACDM model 
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Recently, two issues concerning tiie tiiree-year Wilkinson Microwave Anisotropy Probe (WMAP) 
likeliliood code were pointed out. On large angular scales {I < 30), a sub-optimal likelihood approx- 
imation resulted in a small power excess. On small angular scales {I > 300), over-subtraction of 
unresolved point sources produced a small power deficit. For a minimal six-parameter cosmological 
model, these two effects conspired to decrease the value of Us by ~ 0.7a. In this paper, we study the 
change in preferred parameter ranges for extended cosmological models, including running of ris, 
massive neutrinos, curvature, and the equation of state for dark energy. We also include large-scale 
structure and supernova data in our analysis. We find that the parameter ranges for as, flk and w 
are not much altered by the modified analysis. For massive neutrinos the upper limit on the sum of 
the neutrino masses decreases from AI^, < 1.90eV to < 1.57eV when using the modified WMAP 
code and WMAP data only. We also find that the shift of Us to higher values is quite robust to 
these extensions of the minimal cosmological model. 



I. INTRODUCTION 

The temperature fluctuations in the cosmic microwave 
background (CMB) radiation have proved to be the sin- 
gle most important cosmological observable we have to- 
day, and the high-precision full-sky maps provided by the 
Wilkinson Microwave Anisotropy Probe (WMAP) play a 
very important role in the determination of cosmological 
parameters and preferred cosmological models [1]. One 
of the most important conclusions both from the WMAP 
data alone, and from joint analyses including other cos- 
mological observables, is that a simple six-parameter flat 
ACDM model fits the data very well, and that extended 
models with additional free parameters do not improve 
the fit significantly. 

Because of the great impact of WMAP data, the 3-year 
data analysis from the WMAP team has been subject to 
exhaustive cross-checking. In particular, in refs. [4] and 
3 two noticeable issues with the likelihood code as first 
presented by the WMAP team were pointed out. First, 
the likelihood approximation used between 13 < Z ^ 30 
appears to be inadequate, effectively resulting in a ~ 5% 
power excess in this range compared to an exact treat- 
ment. Second, the amplitude for the unresolved point 
source spectrum used by the WMAP team was found to 
over-subtract the actual contribution in the data, leading 
to a power deficit at high Vs. 



In \§\ , the effect of these discrepancies were studied for 
a minimal six-parameter cosmological model. This was 
done both for WMAP data only, and with additional 
CMB data from small-scale experiments. Their main 
finding was an increase in n^, lowering the significance 
of 71s ^ 1 from ^ 2.7(7 to 2. Oct. 

In this paper we consider the effect on extended cos- 
mological models. We study how the modified WMAP 
likelihood affects the preferred ranges of the running of 
the scalar spectral index, the cosmological neutrino mass 
limits, spatial curvature, and the equation of state for 
dark energy. We have also taken into account large-scale 
structure (LSS) and type la supernovae (SNIa) data sets, 
to see whether the shifts in preferred parameter ranges 
survive a more thorough cosmological analysis. Further, 
we consider whether the shift of to larger values is 
robust to changes in cosmological models and data sets. 

In the next section, we review both the methods and 
data we use. In Section IIIIl we report and comment 
upon our results, before summarizing and concluding in 
Section HVl 
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II. DATA AND METHODS 

A. WMAP data 

The WMAP experiment is a NASA-funded satel- 
lite mission designed to measure the CMB temperature 
anisotropics over the full sky at five frequencies between 
23 and 94 GHz with unprecedented angular resolution 
and sensitivity. These measurements allow for an accu- 
rate determination of the angular CMB power spectrum 
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for angular scales between, say, / = 2 and 800 with three 
years of observations. 

Estimation of this spectrum and the corresponding 
likelihood function is a multi-step process. First, sky 
maps are generated from the raw satellite data, and the 
instrumental noise is estimated. Second, contaminants 
in the form of galactic and extra-galactic foregrounds 
are removed from the sky maps, and severely contami- 
nated regions are removed completely from further anal- 
ysis. Third, the power spectrum is estimated with some 
algorithm, usually trading off computational efficiency 
against accuracy. Fourth, a connection is made between 
the power spectrum and the likelihood. 

These steps are all described in detail for the three- 
year WMAP data release in ref. The main result 
of these efforts is a user-friendly Fortran 90 code that 
for an input power spectrum outputs the corresponding 
likelihood value. In principle, this piece of code may be 
used as a "black box" . 

However, some care is warranted. In particular, two 
points were noted in ref. [2]. First, there is a ~ 5% 
descrepancy between the temperature likelihood approx- 
imation used by the WMAP team and an exact evalua- 
tion for I < 30. Second, there is a ~ 60 /iK^ discrepancy 
between the two spectra observed at 61 and 94 GHz. 
The former is primarily due to estimator approximations 
and secondarily to residual foregrounds. The latter issue 
was later partly explained in terms of an excessive point 
source correction applied to the WMAP spectrum Q . 

In the present paper, we therefore use two versions of 
the WMAP likelihood. The first version is simply the of- 
ficial code as provided on LAMBDA 25]. The second ver- 
sion includes two modifications to this code: At I < 30, 
we replace both the WMAP pixel-based likelihood and 
the pseudo-Cj-based likelihood with an exact Gibbs sam- 
pling based estimator . Then the spectrum amplitude 
of unresolved point sources (relative to 41 GHz) is ad- 
justed from A = 0.017/iK^sr to A = 0.011/xK^sr ^. 
We do not marginalize over the SZ (Sunyaev-Zeldovich) 
amplitude in our analyses. 



B. Other data sets used 

In our analysis we use additional CMB data, data from 
LSS surveys, SNIa data and additional priors on the Hub- 
ble parameter and baryon content of the universe. 



1. Other CMB observations 

To probe a larger range of angular scales in the CMB 
power spectrum we use CMB data from ACBAR [5] and 
BOOMERanG UMM- 



2. Large scale structure 

Large scale structure surveys probe the matter distri- 
bution in the universe by measuring the galaxy-galaxy 
power spectrum Pg(k,z) = {\6g{k, z)\'^) . In the lin- 
ear perturbation regime it is expected that this galaxy- 
galaxy spectrum is proportional to the total matter 
power spectrum, P™, through the simple relation Pg = 
b'^Pm, where b is called the bias parameter. 

There are two galaxy surveys of comparable size, 
namely the 2 degree Field Galaxy Redshift Survey (2dF) 
and the Sloan Digital Sky Survey (SDSS) [l0|. In our 
analysis we use data from both these surveys. 



3. Type la supernovae 

Probing the luminosity-redshift relation of SNIa is one 
of the most direct measurements of cosmological expan- 
sion, and thus one of the most powerful pieces of evidence 
for the existence of dark energy. In our analysis we use 
SNIa data from the Supernova Legacy Survey (SNLS) 
which is a dedicated SNIa survey currently includ- 
ing 71 SNIa in the redshift range z = 0.2 — 1. 



4. Additional priors 

In addition to the CMB, LSS and SNIa data sets men- 
tioned above, we impose priors on the Hubble parameter, 
the baryon content in the universe, and the position of 
the LSS baryonic peak. 

From the Hubble Space Telescope Key Project (HST) 
we have ado pte d a prior on the Hubble parameter of h = 
0.72 ± 0.08 [141. The constraint on the baryon density 
today was chosen to be flbh"^ = 0.022 ± 0.002 from Big 
Bang nucleosynthesis (BBN) [l3t[ll|. 

From the detection of baryonic acoustic oscillations 
(BAO) in the sample of luminous red galaxies (LRG) in 
the SDSS survey [l5| it is also possible to put a constraint 
on the combination of parameters 



A 



BAO 



H{z) 



1/3 



(1) 



where Dm{z) is the comoving angular diameter distance. 
The BAO constraint can then be written as 



A 



BAO 



0.469 



f— ) 
Vo.98/ 



-0.35 



(1 + 0.94/,) ± 0.017, (2) 



where the fit to the neutrino fraction, ~ ri,/r2„i is 
given by [IB]. For z we adopt the redshift of a typical 
LRG in the SDSS sample, z = 0.35. For models with 
non-zero as , we substitute in eq. ^ with an effective 
Us given by 



gff(fci) 



dlnP 



c?ln fc 



+ 1 = ns{ko) + as ln(fci/fco). 



k—ki 



(3) 
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where P is the primordial power spectrum given by 

InP = In^ + {ns - 1) Hk/ko) + as/2ln{k/kaf , (4) 

and ko is set to O.OSMpc"^ We use fci = 0.01/iMpc~\ 
which is approximately the scale associated with the 
baryonic peak in the LRG sample. Our results are ro- 
bust to changes of fci around this value, as the preferred 
values of as are small in the models considered here. 



C. Parameter estimation 

The parameter estimation process is based on the pub- 
licly available CosmoMC code using the data and 
likelihoods described above. 

For each model, we compute the corresponding pa- 
rameter confidence intervals using three different com- 
binations of data sets, named A, B and C (see Table 
Data set A includes the three-year WMAP data only; 
data set B includes also CMB data from ACBAR and 
BOOMERanG; and data set C also LSS and SNIa data 
sets and priors from HST, BBN and BAG. 

As a basic six-parameter cosmological model we use the 
parameters {fif,/i^, fim, log(10^*'^s), /i, n^, r}. The exact 
parameter definitions are given by the CosmoMC code, 
fift is the ratio of baryons to the total energy density; h is 
the Hubble parameter today; flm is the ratio of matter to 
the total energy density today; Ag sets the amplitude of 
the primordial fluctuations; Ug is the tilt of the primordial 
power spectrum; and r is the reionization optical depth. 

We then extend the six-parameter model by adding 
the parameters as,r, M^,nk,w one by one (except for 
as and r which are added simultaneously). Here, as is 
the running of n^, defined in eq.(l4]); r is the ratio of ten- 
sor to scalar fluctuations; Mi, is the sum of the neutrino 
mass eigenstates, M^, — 'Y^rrii, = 93.14i7iy/i^; fi^ is the 
amount of spatial curvature; and w is the equation of 
state parameter for dark energy. Finally, we vary all 11 
parameters simultaneously. 

For all combinations of data and parameter sets, we 
carry out a similar analysis both with the standard 
WMAP likelihood code as provided and with the two 
modifications described above. 



Data set 


Observations included 


A 


WMAP 


B 


WMAP + ACBAR + BOOMERanG 


C 


WMAP + ACBAR + BOOMERanG + SDSS 




+ 2dF + SNLS + HST + BBN + BAO 



TABLE I: The different combinations of data sets used in this 
analysis. 



III. RESULTS 
A. Minimal six-parameter model 

We start with the simple six-parameter model hav- 
ing the free parameters {^ih/i^, rim, log(lO^'^As), h, ng, r}. 
This was first done in ref. Q for the combinations of data 
sets A and B (see TableHl. For Us they reported a ^ O.Tcr 
shift to higher values when applying the modified anal- 
ysis. The other parameters were only subject to small 
shifts of their mean values. 

We repeat this analysis here, but also include data set 
C in the analysis. Consistent with ref. [3| we find that 
only Us is notably affected by the modified likelihood, 
and the shift of Us also remains for data set C, in which 
case the value for rig changes from Us = 0.961 ± 0.014 to 
ns = 0.971 ±0.014. This corresponds to a shift of ^ 0.8cr, 
and weakens the significance of 7^ 1 from ~ 2.9 to 
^ 2.1 for this data set. As can be seen from refs. 0,0], 
the low-Z and point-source corrections contribute almost 
equally much to the shift in n^. In ref. Q they found a 
mean value of Us of 0.961 when applying only the low-/ 
corrections. 

That the shift of ris survives when adding LSS data 
is not very surprising; is less sensitive to LSS than 
to CMB data because of the larger dynamic range and 
higher precision of the latter observations. The resulting 
values of are summarized in Table [Hi 



Data set WMAP code 


Modified code 


A 


0.954 ± 0.016 


0.966 ±0.016 


B 


0.958 ±0.016 


0.969 ±0.016 


C 


0.961 ±0.014 


0.971 ±0.014 



TABLE 11: Results for ris in a six-parameter model. The val- 
ues in the second column are found using the WMAP likeli- 
hood code, while the values in the third column are calculated 
using the modifications described in section III Al All errors 
are la. 



B. Running of spectral index 

Next we consider how the modified WMAP likelihood 
affects the constraints on ag. The simplest inflationary 
models predict an as that is slightly different from 0, 
and thus information on as can provide us with valu- 
able information on inflationary mechanisms. Following 
Spergel et al. [l[, we marginalize over the ratio of tensor 
to scalar fluctuations, r, since models with negative as 
often correspond to large tensor modes. 

Our results are summarized in Table [1111 We find that 
the likelihood modifications have no major effect on the 
constraints on as, but we observe a small increase of 
~ 0.2cr in the significance of as 7^ 0. 
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Data set 


WMAP code 


Modified code 


A 


-0.050 ±0.027 


-0.052 ± 0.027 


B 


-0.052 ±0.026 


-0.056 ±0.025 


C 


-0.013 ±0.020 


-0.014 ±0.019 



TABLE III: Estimated values for from the WMAP like- 
lihood code and our modified code. The other parameters, 
including r, are marginalized over. 



C. Massive neutrinos 

Another natural extension of the minimal six- 
parameter model is the addition of massive neutrinos. 
This is motivated by observations of neutrino oscillations, 
which show that neutrinos indeed are massive. 

Because of their low mass, neutrinos act like a warm 
dark matter component in the universe. Given the energy 
fraction of massive neutrinos today, U,^, one can easily 
find a limit on the sum of the neutrino masses, M^, by 
the relation Vt^K^ = M^/93.14eV. (See ref. [l| or [3 
for a review of the cosmological properties of massive 
neutrinos.) 

At present, the best upper limits on the absolute mass 
scale of neutrinos come from cosmology. The current 
cosmological 95% C.L. limits range from < 0.17eV 
[20| . relying on extensive use of different data sets, to 
< 2.0eV for WMAP data only [H, EH, H^. In ref. 
[2H they also pointed out that it will be difficult to push 
the upper limit on r^i^/i^ much below fi^/i^ ^ 0.017 using 
CMB data only. This corresponds to a neutrino mass 
limit Ml, < 1.5eV. For smaller neutrino masses, neutrinos 
will still be relativistic at the time of recombination, and 
thus the effects of the neutrino masses will not be fully 
revealed in the CMB power spectrum. 

In our analysis we assume three species of massive neu- 
trinos with degenerate masses. The assumption of degen- 
erate masses has been shown to be verygood for the mass 
regime that we are working in here [23j. The resulting 
neutrino mass limits are summarized in Table IIVI We 
see that when using WMAP data alone, the upper limit 
on Ml, is significantly improved by the modified analy- 
sis, and that we are approaching the limit of how tight 
constraints on M^, we can find from CMB data alone. 
If we analyze the data applying only the low-Z correc- 
tions to the WMAP code, the neutrino mass limit be- 
comes Ml, < 1.69el^, which shows that both the low-^ 
and point-source corrections are important also for this 
model. 

Our improved Mi, limit can be understood by the slight 
degeneracy between the Mi, and Hs parameters, in that 
a larger value of Us provides less space for a large Mi,. 
This can be seen from the contour plot in the Mj^-rig 
plane in Figure [1] The degeneracy can be understood 
by the fact that both rig and Mi, have impact on the 
small-scale behavior of the CMB power spectrum. 

By definition, Us sets the tilt of the primordial spec- 
trum. If Ml, is of order ^ 2eV, it will affect the power 



on scales smaller than / ~ 300. This happens because 
the perturbations of the gravitational potential on scales 
smaller than this are suppressed by neutrino free stream- 
ing, which in turn boosts the acoustic oscillations ^231. 
As Ml, increases and more of the dark matter consists 
of massive neutrinos, this boost of small scale power also 
increases. Therefore, a large value of increases the 
power on small scales, leaving less room for to add 
further power without coming in conflict with data. Mi, 
also affects the heights of the peaks on larger scales, but 
this can to a large extent be compensated for by adjust- 
ing the values of flm and fl^h^. 

It is interesting to notice that the upper limit on Mi, 
actually weakens when we add small scale CMB data 
sets. From Table HIl we also see that the preferred value 
of 71s also increases when these data sets are included. 
This indicates that the small-scale power from ACBAR 
and BOOMERanG is higher than what one would expect 
from the WMAP data. When adding massive neutrinos 
to the minimal six-parameter model this increment in 
small-scale power can be partly accommodated by in- 
creasing Ml, instead of Ug . This will increase the small- 
scale power without altering the fit to the large scale 
spectrum. 

We also see that when we add LSS and SNla data, the 
neutrino mass limit is no longer affected by the modified 
WMAP analysis, as the additional constraints on Mi, are 
mainly determined by LSS data. In the LSS power spec- 
trum neutrino free-streaming will suppress small scale 
power, and a larger will in this case allow for a larger 
M,y. This effect then cancels out the improved Mi, limit 
found from the WMAP data. 



Data set WMAP code Modified code 



A 


1.90eV 


1.57eV 


B 


2.13eV 


L72eV 


C 


0.45eV 


0.45eV 



TABLE IV: Estimated 95% C.L. upper hmits on M„. 



D. Spatial curvature 

Next we add spatial curvature, Jlfe, to our six- 
parameter model. However, altering the geometry of uni- 
verse mainly affects the positions of the CMB acoustic 
peaks, while the likelihood modifications mostly concern 
the amplitude and tilt of the power spectrum. A priori, 
one would therefore not expect any significant changes in 
fifc. And as seen in Table IV| this is indeed the case. For 
all data sets, f2fe = is within ^ Ict, in good agreement 
with the results from the WMAP team [l[ . 

The large improvement of the limits on fJfc for data set 
C can to a large extent be understood by the well-known 
degeneracy between $7^ and /i, where negative values of 
rife can be accommodated by a small h. Therefore, when 
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0.2 0.4 0.6 0. 



FIG. 1: 68% and 95% C.L. contours in the Mi,-ns plane for 
the modified likelihood code and a seven-parameter model 
with free neutrino mass, using WMAP data only. There is 
a slight degeneracy between the two parameters, and larger 
values for Ua will put tighter upper limits on M^. 



imposing the HST prior on h, the allowed range of f2fe is 
significantly constrained. 

Data Set WMAP code Modified code 



A 
B 
C 



U.UO( _o.056 

_n n5fi +0 052 

U.UOD _o.062 



U.UO/ -0.057 



-0.055 



+0.048 
-0.055 



-0.005 ±0.007 -0.006 ±0.007 



TABLE V: Estimated values for Q.k 



E. Dark energy equation of state 

The nature of dark energy is one of the major questions 
in cosmology today. In the minimal six-parameter model, 
the dark energy is assumed to be a cosmological constant 
with w = — 1, and this has been shown to agree well with 
current cosmological data (see, e.g., ref. Here we 

test whether the modified WMAP likelihood code alters 
the preferred values for w. In the following analysis we 
assume that w is independent of redshift. 

The main effect of w on the CMB power spectrum is 
to shift the position of the acoustic peaks by altering the 
expansion history of the universe. Therefore we would 
not expect the limits on w to be much affected by the 
new WMAP likelihood analysis. Still, from Table IVll we 
notice a small shift of order ~ 0.2(t to smaller values of 
w when using CMB data only. This happens because a 
smaller w will enhance the late integrated Sachs- Wolfe ef- 
fect, which results in a suppression of large-scale fluctua- 
tions in the observed CMB power spectrum. The slightly 
smaller preferred value of w in the modified analysis is 
accompanied by small changes also in h and to shift 



the peaks back in position. As more data is added, we see 
that the modified analysis has no effect on w anymore, 
and that a cosmological constant still fits the data. 



Data set WMAP code Modified code 

A -0.98 ±0.41 -1.05 ±0.39 

B -0.97 ±0.41 -1.03 ±0.37 

C -1.00 ±0.07 -1.00 ±0.07 

TABLE VI: Estimated values for w. 



F. 11-parameter model 



Finally we vary 



all 



11 



parameters 



{Vtbh'^, Q.m-, \og{lQ'^°As),h, Us, r, a^, r, M^, ilk, w} si- 
multaneously. The results for the parameters Hs, as, 
ilk, w and Ml, are shown in Table rvTTl Here we see that 
all effects found in the more restricted models above are 
also present in this extended model. 

For Hs the likelihood corrections still result in a mean 
value that is ~ 0.02 larger than with the original like- 
lihood. However, as the uncertainty in ris is increased 
(mainly due to a degeneracy with a^), this shift is not 
as statistically significant as it was in the six-parameter 
model in subsection IIII Al Further, for data set C we 
see that the preferred value of Ug is in fact not much 
changed by the modified analysis. Rather, the power 
spectrum changes are accommodated by a slightly lower 
value of as to account for the smaller power for low Vs 
in the CMB power spectrum. 

For as , there are no significant changes of the preferred 
values by the 11-parameter model, and we see that ag = 
is still consistent with the data. 

The same is the case for flk- Here the modified anal- 
ysis shifts the preferred values to slightly more negative 
values, but with all data sets, fifc = remains well within 
la of its mean value. 

Also for Ml, we find the same effects as above. When 
using CMB data only, the upper limit on M^, is improved 
by the new WMAP likelihood analysis. As expected, the 
limit becomes weaker for the extended parameter set, but 
from WMAP data alone we still find a M,, limit that is 
better than what is found in earlier papers, even with 
this large parameter space. Also, we see that the Mi, 
limit still becomes weaker by adding the ACBAR and 
BOOMERanG data sets. 

For the dark energy equation of state, we still find 
that the modified analysis prefer slightly lower values for 
w. In this extended model the preferred values of w are 
shifted to lower values than in the more restricted model 
in subsection llll El which is mainly accommodated by the 
preferred negative value of flk (as w and Qk shifts the 
positions of the acoustic peaks in opposite directions). 
Still w = —1 remains within la for all of the data sets. 

In Figure [21 we show the best-fit power spectrum for 
the 11-parameter the model both from the analysis with 
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the standard WMAP likelihood code and the new point- 
source and low-^ corrected likelihood analysis. We see 
that the discrepancy is most notable for Z < 100. 



Parameter 


WMAP code 


Modified code 


Data set A 


Ms (0.05) 


0.863 ±0.047 


0.880 ±0.046 


Ols 

w 


-0.051 ±0.029 
-0 019 

U.Uia _o.053 

-1.43 ±1.09 


-0.050 ±0.029 

-0 027 +°-°'^^ 
-1.53 ±1.12 




<2.09eV @ 95% C.L. 


<1.66eV @ 95% C.L. 


Data set B 


ns{0.05) 


0.859 ±0.042 


0.875 ±0.041 


Cts 

w 


-0.055 ±0.027 
-0 010 +""*'^ 

U.UiU _o.050 

-1.33 ±1.02 


-0.055 ±0.027 

n ni« +0.048 
-0.018 „m,53 

-1.43 ±1.10 




<2.33eV @ 95% C.L. 


<2.02eV @ 95% C.L. 


Data set C 


ns{om) 


0.954 ±0.038 


0.954 ±0.036 


Cts 


-0.003 ±0.028 


-0.012 ±0.026 


w 


-0.001 ±0.012 
-1.05 ±0.09 


-0 001 +''°" 

U.UUl „Q 013 

-1.05 ±0.09 




<0.51 eV @ 95% C.L. 


<0.52eV @ 95% C.L. 



TABLE VII: Parameter results for the model that includes 
free Us, r, M^, w and Q^. 



IV. CONCLUSION 

In refs. 

and [S*], two modifications to the three-year 
WMAP likelihood were presented. Relative to the power 
spectrum presented by Hinshaw [1], they found a power 
deficit for low values of I due to inaccurate likelihood 
approximation, and a small power excess for high values 
of I due to over-subtracted unresolved point sources. 

The impact on the inferred cosmological parameter in- 
tervals from these corrections in a minimal six-parameter 
model was studied in ref. [1| using CMB data only. Their 
single most important result was an increase in the pre- 
ferred value of rig, lowering the significance of ris ^ — 1 
from ^ 2.7a to 2.0a. In the present paper, we have 
extended that analysis to also account for cosmological 
models including a non-zero running of Hs, massive neu- 
trinos, curvature and w —1. We have also added LSS 
and SNIa data sets to our analysis to see if the parameter 
shifts induced by the modified WMAP analysis survive 
when adding more data sets. 

We found that the shift of Ug to larger values survives 
when we add LSS and SNIa data. However, when we 
apply all data sets in the full 11-parameter model, Us 
is not affected much by the modified WMAP analysis 
anymore. This is mainly due to the extra freedom with 
as. 
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FIG. 2: Upper panel: The CMB power spectrum for the 
best-fit 11-parameter model with the WMAP likelihood code 
(solid blue line) and the low-Z and point source corrected like- 
lihood code (dashed red line). The gray shading shows the 
cosmic variance around the blue line. Note the small dis- 
crepancies at both low and high Vs. Lower panel: The rela- 
tive difference between the two power spectra, (^modified _ 
^WMAP^^^WMAP 



For the extended models, we found that the preferred 
values of a^, and w are not significantly affected by 
the modified analysis. When including massive neutrinos 
we found that the upper limit on when using WMAP 
data alone was reduced from AI^, < 1.90eV to M^, < 
1.57eV. A similar improvement in the limit could not 
be observed when adding LSS and SNIa data, since the 
higher preferred value of Us will allow for larger neutrino 
masses in the LSS power spectrum. 

Since the initial publication of the two re-analysis 
papers, refs. and @, and the present paper, the 
WMAP team has released a new version of their like- 
lihood codej26j that implements the suggested \ow-l cor- 
rection and a revised point-source correction. Using this 
updated likelihood code, we find Ug = 0.959 ± 0.016 
for the six-parameter model and the WMAP data only. 
Including massive neutrinos, this code gives an upper 
bound of < 1.75eV from WMAP data only The 
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difference is due to the point source amplitude and cor- 
responding error adopted by the WMAP team, which do 
not match perfectly that of ref . 3 ■ Unfortunately, full 
details on the WMAP approach are not currently avail- 
able, and final assessment of this issue must therefore 
await the release of the revised WMAP3 papers. 

To conclude, the modified analysis does not strengthen 
the case for non-standard cosmological parameters, and 
the standard flat ACDM model still provides an excellent 
fit to data. 
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